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1 Emptiness formation probability at large distance 



The Hamiltonian of the XX Z spin-| Heisenberg chain is given by 

M 

H=Y1 + «+i + ~ !)) • (1-1) 

m=l 

Here A is the anisotropy parameter, am V ' z denote the usual Pauli matrices acting on the quan- 
tum space at site m of the chain. The emptiness formation probability r(m) (the probability 
to find in the ground state a ferromagnetic string of length m) is defined as the following 
expectation value 

T{m) = m\l^±\i> 9 ), (1.2) 

k=l 

where \ifj g ) denotes the normalized ground state. In the thermodynamic limit (M — > oo), this 
quantity can be expressed as a multiple integral with m integrations |l], |2[ ||, ||, ||] . Recently, in 
the article Q, a new multiple integral representation for r(m) was obtained. It leads in a direct 
way to the known answer at the free fermion point A = ||, in particular using a saddle point 
method, and to its first exact determination outside the free fermion point, namely at A = \ 



1 10]. 



The purpose of this letter is to present the evaluation of the asymptotic behavior of r(m) 
at large distance m, in the massless regime — 1 < A < 1, via the saddle point method. We find 

logr(m) 7T 1 f duj sinh %(ir - Q cosh 2 %f 

lim & I ' = log - + - / ^ ^- (1.3) 

m-*oo m z Q 2 J uj sinh W- sinh ^ cosh ujC 

R-iO 2 2 

where cos^ = A, < £ < it. If £ is commensurate with ir (in other words if is a root of 
unity), then the integral in (|1.3j ) can be taken explicitly in terms of tp- function (logarithmic 
derivative of T-function). In particular for Q = ^ and C = f (respectively A = and A = 1/2) 



we obtain from (1.3) 



log r(m) 1 , 
lim \ ' = -log 2, A = 0, 

(1.4) 

log r(m) 3 1 
hm — = -log3-31og2, A = -, 

m-*oo m z 2 2 

which coincides with the known results obtained respectively in [0, |8|, § and in [11, 1C]. For the 



particular case of the XXX chain (A = 1, £ = 0) the asymptotic behavior can be evaluated 
also by the saddle point method and it is given by 

„ m !^q™> _ log ( T JIM>) K log ( ^5991), (1.5) 

m-^oo m z \ T(|) I 
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which is in good agreement with the known numerical result log(0.598), obtained in |12|| , 

Below, we explain the main features of our method. A more detailed account of the proofs 
and techniques involved will be published later. 



2 The saddle point method 

The multiple integral representation for r(m) obtained in H can be written in the form 



'W = (afSE?)" '(?) 2 J^ F (^ 

V 

" sinh | (A Q - A,) - / si nh(A a - |) sinh(A a + f> . 

W sinh(A a - X b - i() sinh( A a - A 6 + <) ±± ^ cosh f A a J A • J 



with 

/ r \ ~i \ , 1 / — zsinC \ 

F({A},m = hm — det m — — — . - . ; — . 2.2) 

Ci,.^— § nsinhfe-eft) Vsmh(A i -4)smh(A,-e fc -<)y 1 j 

Here the integration domain D is — oo < Ai < A2 < • • • < X m < 00. 

Following the standard arguments of the saddle point method we estimate the integral (2.1) 
by the maximal value of the integrand. Let {A'} be the set of parameters corresponding to 
this maximum. They satisfy the saddle point equations and for large m we assume that their 
distribution can be described by a density function p(X'): 

piX'j) = hm -—i — . (2.3) 

Thus for large m, one can replace sums over the set {A'} by integrals. Namely, if /(A) is 
integrable on the real axis, then 



71 j=i 

m — > 00. (2.4) 



m 3=1 'V /v fc J -°° /x n h 



Due to (2.4) it is easy to see that in the point A' l5 . . . , X' m the products in the second line of 
(2.1) behave as exp(c m 2 ). 
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Our goal is now to estimate the behavior of the term F({X'},m). To do this we factorize 
the determinant in (|2.2|) as follows for large m: 



det r 



-i sin £ 



with 



sinh(A^. - sinh(Aj- - & - <) y 



ftr m j Sil1 2 C f OR , 

smn(A — z^J) smn(A + zi, J 



Indeed, for m — > oo one has 

det m I djfc - - — tttt" I det 



( 



2mmp{\' k )J m V2Csinh|(A;. -e 



det r 



jm/o(Aj) 2Csinh^( 



6 m v2Csinh^(A;. -^ fc ) ^ 2vimp{\\) 2Csinhf (A{ - 1 ' j 



2Csinhf(A;.-6) i-oc 27ri 2Csinhf(/i-e fe ; 



1 \ m / sin t 

— I det m ' 



2tt ; I sinh(A' - &) sinh(A' - & - <) 



Here we have used the fact that the function i/2Q sinh ^(Aj — £) solves the Lieb integral equation 
for the density of the ground state of the XXZ magnet [jl3| (and we have used the notations 
of §) . The second determinant in the r.h.s. of ( |2.5| ) is a Cauchy determinant, hence, 

m 



II cosh m ?A' a 

a=l 



The behavior of the determinant in fl2,8|) can be estimated via Hadamard inequality 

|det m (o 3 -fc)| < (max\a jk \) m m^ . (2.9) 

applied to the above determinant and to the determinant of the inverse matrix , which shows 
that 

1 / K{X' j -X' k )\ 

lim — logdet m 5 jk - \ J ../Y = 0. 2.10 
m-^oo m z \ ImmpyX, ) J 
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The last equation means that det m ( Sjk — K(\ ! - — X' k ) /2mmp(X' k ) j does not contribute to the 
leading term of the asymptotics. Hence, it can be excluded from our considerations. 

Thus, up to subleading corrections of the exponential type the emptiness formation proba- 
bility behaves as 



r(m) 



, m 

c; 



m — > oo, 



(2.11) 



with 



S = S({X'}) = — J> 



sinh 2 f(A^-A^ 



, v sinh(A' a -A^-iC)sinh(^-A' b + iC) > 
+1 y ^ ( sinh ( A a " l 'C/2) sinh(A' a + iC/2) \, 



a=l 



cosh 2 |A' a 



(2.12) 



Here the parameters {A'} are the solutions of the saddle point equations 

dS 



dX' : 



0. 



(2.13) 



In our case the system ( 2.13 ) has the form 



2tt 7tA'- 

— tanh — ^ - coth(A' - i(/2) - coth(A'- + i(/2) 

c c 



1 m /2vr 7T \ 
- E (— «>th -(A$ " A' fc ) - coth(A< - A' fc - *C) - coth(A;- - X' k + . (2.14) 



Using ( |2.4[ ) we transform ( 2.14| ) into the integral equation for the density p(X) 



^ tanh ^ - coth(A - i(/2) - coth(A + i(/2) 

= V.P. J ( — coth - (A - At) - coth(A - p-iQ- coth(A - /i + <) ) /?(/i) d/i. (2.15) 

Respectively the action So takes the form 

' sinh( A - i(/2) sinh(A + i(/2) \ 



SV, 



dXp(X) log 



cosh 2 ?A 



+- 



dpdXp(X)p(p) log 



sinh 2 f(A-/u) 



sinh(A — p — iQ) sinh(A — p + z£) 



(2.16) 
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Since the kernel of the integral operator in (2.15) depends on the difference of the arguments, 
this equation can be solved via Fourier transform. Then 



p(u) 



Making the inverse Fourier transform we find 

P(A) 



p(X) dX 



coshf^ 



cosh 
cosh ujC 



(2.17) 



C\/2cosh^' 



(2.18) 



which obviously satisfies the needed normalisation condition for density (integral on the real 
axis equals one). It remains to substitute (2.17), (2.18) into ( 2.16p , and after straightforward 
calculations we arrive at 



So - g 



R-iO 



duj sinhf (tt - C)cosh 2 s£ 
w sinh^sinh^cosh^C 



(2.19) 



Thus, we have obtained (1.3). 

In the case of the XXX chain (A = 1) one should rescale Xj — > (Xj, £j — * ££j in the original 
multiple integral representation ( p.l[ ) for r(m) and then proceed to the limit £ — ► . The 
remaining computations are then very similar to the ones described above, therefore we present 
here only the main results. The behavior of r(m) is now given by 



r(m) 



m 



oo. 



(2.20) 



The action Sq in the saddle point has the form 

'(A-i/2)(A + «72) 



So 



log 



1 

+ 2 



cosh 2 7rA 



dfid\p(X)p(fi) log 



p{\) dX 

sinh 2 7r(A — p) 



(A — p, — i)(X - p + i) 
The analog of the integral equation ( |2.15| ) in the XXX case is 



(2.21) 



2ir tanhvrA 



2A 



A 2 + 



V.P. 



2tt coth 7r(A — p) 



2(A - p) 
{X-pf + l 



p(p) dp. 



The solution of this equation is 



P(A) 



cosh ^ 
\/2cosh7rA 



(2.22) 



(2.23) 



Substituting (p.23|) into (2.21) we finally arrive at (|1.5|). 
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